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A detailed investigation of the temperature dependence of the spatial string tension ds 
in SU[2) gauge theory is presented. A sustained performance of 3 GFLOPS on a 64K 
Connection Machine CM-2 equivalent has been achieved. Scaling of ds between f3 = 
2.5115 and f3 = 2.74, on large lattices, is demonstrated. Below the critical temperature, 
Tc, (Ts remains constant. For temperatures larger than 2Tc the temperature dependence 
can be parametrized by <Js{T) = (0. 369 it 0.01 4) (T)T^ , where g{T) is a 2-Ioop running 
coupling constant with the scale parameter determined as A^^ = (0.076 it 0.013)Tc. 

Keywords: Lattice gauge theory. Finite temperature QCD, Dimensional reduction. 
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1. Outline 

After a brief introduction into lattice gauge theories where we concentrate on the 
aspects of computational relevance, we discuss our implementation of the updating 
algorithm on the Connection Machine CM-2. This is followed by a presentation of 
the physical background, motivation, and a summary of related results, obtained by 
previous authors (Section 3). In Section 4 technical aspects of our data evaluation 
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procedure are discussed. The paper concludes with physical results (Section 5). 



2. Simulation of SU(2) Gauge Theory 
2.1. The ''mathematical" problem 

The lattice regularization is the only known non-perturbative method to treat quan- 
tum field theories quantitatively. A space-time lattice with lattice spacing a and 
volume (N^ x Nr)a'^ is used. Physical observables can be extracted in the contin- 
uum limit (a 0). An element of the gauge group, responsible for interactions, is 
mapped to each link connecting two lattice sites. In the physical theory of inter- 
est, the theory of strong interactions, QCD, this group is SU(3). Because of the 
similarity of SU(3) vacuum structure with the one of SU(2) we simulate SU(2), in 
order to minimize the computational task. 

A set of group elements (link variables) C = {Ufi(n)} with //G{l,...,4},n = 
(rii, n2, "3, "4), rii £ {1, . . . , N„}, 714 £ {1, . . . , Nr}, and U^{n) £ SU (2) is called 
a configuration. Each configuration represents a possible field distribution on the 
finite lattice. In the path integral formulation of quantum field theories, expectation 
values of operators (O) (observables) can be calculated as an average of the values 
O takes on all possible configurations, weighted by a Boltzmann measure: 



with l3S being a discretized version of the continuum action of the Euclidean theory. 

The simplest possible gauge invariant choice with the desired continuum limit 
is the so-called Wilson action: 



where Ufi^(n) denotes the ordered product of all link variables around the elemen- 
tary square (plaquette) emanating from site n into the positive oriented //i^-plane. 
In this choice /3 = 4/(jf^ with g being the bare (lattice) coupling parameter. By 
changing /3 one can control the lattice spacing a. Since nonabelian gauge theories 
are asymptotically free, the continuum limit corresponds to infinite /3. The dimen- 
sionful quantity a can me measured in units of some physical scale, e.g. a mass like 
the square rott of the string tension ^/a, or the critical temperature Tg. 

The (high dimensional) integral Eq. 1 can be calculated numerically by use 
of stochastic (Monte-Carlo) methods. As most possible configurations are strongly 
suppressed by their Boltzmann weight factors, one better concentrates on the es- 
sential configurations. This is achieved by the technique of importance sampling: 
Configurations are generated according to the probability density 
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After generation of a suitable number of such (independent) configurations tlie 
observable of interest is calculated as an average of the values of the corresponding 
operator on the ensemble of generated gauge fields. 

The ensemble is produced by a Markov process: Each configuration only de- 
pends on a finite number of preceding configurations: €„ = /(Cn-i , C„_2, • • • , Cn-m) 
with m being the order of the Markov-chain. A criterion which guarantees the right 
equilibrium distribution for the Monte-Carlo simulation is the so-called detailed 
balance condition 

dP{Cn-l)W{Cn-l ^ Cn) = dP(C„)W(C„ ^ (4) 

with W(Cn-i Cn) being the appropriate transition probablility. In addition to 
this requirement the algorithm has to be ergodic, i.e. there has to be a nonvanishing 
probability for each possible configuration to be generated within a finite number 
of steps. For a limited simulation time this means that the trajectories should cover 
the phase space as uniform as possible. 

2.2. Implementation on the Connection Machine CM-2 

There is some freedom in choosing W(Cn-i €„) according to Eq. 4. In our 
implementation we use a combination of two algorithms for updating the gauge 
fields. One of them is the heatbath algorithm, which allows a change in the action 
with transition probability 

W(C„.i^C„)(xeM-f^S(C„)) . (5) 

The necessary group measure is generated by the method proposed by Kennedy 
and Pendleton,^ which is suitable for SIMD architectures. This method is combined 
with the deterministic microcanonical Creutz-overrelaxation^ that provides a faster 
movement of gauge fields in order to fight critical slowing down. Both algorithms 
are mixed stochastically to a hybrid form with the probability of heatbath sweeps 
being varied from 1/4 to 1/14, depending on /3, Nr, and our experience. 

Due to the locality of the Wilson action (Eq. 2), the implementation on a serial 
machine is straight forward: One changes all AN^Nr links one at a time, either in 
a fixed order or randomly. After each link has been updated once (in the average) 
one has carried out a sweep. On a parallel machine one has to be more careful since 
a whole set of links is updated at once. Because of the Markov condition links that 
are updated at the same time must not depend on each other. With the Wilson 
action this means that two links being updated at once must not share a common 
plaquette. The simplest way to provide this is a red/black division of the lattice. 
Lattice sites with even sum of coordinates are called red, the others black. One 
loops over red/black, and the direction All of these eight subsets include no links 
sharing a common plaquette. 

We implement the algorithm on a data parallel programming model (CMFor- 
tran) which is the only possibility on SIMD architectures but also useful for many 



implementations on MIMD machines, such as the CM-5. Each CM-2 has some 
fixed number of physical processors Npp. In order to allow portability of programs 
between Connection Machines with different Npp, and to support scalability of 
codes, physical processors are hidden to the Fortran programmer behind the (soft- 
ware) concept of iJirtaa/ processors (VPs): Each site of a parallel array corresponds 
to a VP. For the programmer all (parallel) commands are performed on all VPs at 
one and the same time. Internally, a certain number of VPs is mapped to one PP 
(therefore, Nyp > Npp to avoid latency of parts of the machine), and the looping 
over these VPs is automatically done by the compiler. For instance, VPs can be 
thought of as sites in a space-time lattice. Then, the canonical data structure for the 
configuration {Ufi(n)} would be: JJ(4:,Ncr,Ncr,Ncr,Nr), where U is an SU(2) matrix, 
with the first dimension being serial (denoting the direction /j, of the link), and the 
other dimensions — being parallelized over — addressing the lattice sites n. 

Implementation of the red/black mask means that the computations are only 
performed on half of the lattice volume at once. Thus, half of the machine would 
be inactive. In order to establish load balancing it is advantageous to declare 
V(4:,2,Ncr,Ncr,Ncr,Nr/2), instead. The second index is serial and labels red/black 
sites. This reduces the number of virtual processors to half of the lattice volume 
and guarantees that all processors are active all the time. 

On the Connection Machine, there are still two problems with this kind of data 
layout. One is excessive memory requirement: Although they are not needed all the 
time (VP looping), all temporary arrays (as they are needed for matrix multiplica- 
tions etc.) have parallel dimension equal to Nyp, which is half the lattice volume. 
So each temporary array consumes one eighth of the space of the configuration itself. 
The second problem is communication: Communication between VPs belonging to 
the same PP is much faster than communication between different PPs. The overall 
speed of a communication step (like sending a link to the left neighbor) is deter- 
mined by the speed of the slowest communication channel involved. In CMFortran 
for CM-2 the programmer has only little infiuence on the mapping of VPs to PPs. 
So a mixture between on- and off-chip communication in one step can easily occur, 
depending on lattice, and machine geometries. From this point of view it is wise to 
emphasize on-chip communication by adopting the machine size in such a way to 
the computational problem that the VP-ratio {Nyp / Npp) is large. 

One possibility to circumvent these problems is blocking the lattice into hyper- 
cubes (of 2-* lattice sites each) by hand: U(4,16,#<,/2,#<,/2,#<,/2,#^/2) The second 
index labels the address within each hypercube (lowest bit: first direction etc.). The 
serial loop is now over 64 instead of eight subsets. Thus, memory requirement for 
temporaries is reduced by a factor eight! For one half of all communications a high 
speed is guaranteed since the sites are located on the same VP and "communica- 
tion" is programmed by hand. The disadvantages are higher complexity of Fortran 
code, limitation to lattices of even extents N^, and Nj, and a lower VP-ratio (by 
this factor eight). 

For an 8K CM-2 (s/zcewzse model) like the one in Wuppertal Nyp > Npp = 2"'^'^. 



This translates into 2^^ lattice sites. Saturation in the sense that the GFlops rate 
remains constant within 10% is reached for 2^'^ lattice sites (e.g. on a 32^ x 2 lattice). 
On the 16K CM-2 (fieldwi.se model) at HLRZ Npp = 2^^. For saturation the lattice 
must have at least 32'* sites. Using the above kind of algorithm we reach a sustained 
all-over speed (updating and measurement) of 1.5 GFlops/64K {fieldwtse model) at 
the HLRZ or 3 GFlops/64K [sUcewise model) at the Wuppertal and SCRI CM-2s. 
The latter rate has to be compared with 4.5 GFlops/64K in the SU{?>) case where 
the computation/communication ratio is larger. To reach the above performance 
the matrix multiplications have been programmed in assembler language. 

3. Physical Background 

An understanding of nonabelian gauge theories is important for the physics of the 
quark-gluon phase of matter to be studied in heavy ion collisions. The temper- 
ature T corresponds to the inverse temporal extent of the lattice: aT = l/Nj. 
Non-abelian SU(N) gauge theories in (3-|-l)-dimensions are known to undergo a 
deconfining phase transition at high temperature. The physical string tension, 
characterizing the linear rise of the potential between static quark sources with 
distance, decreases with increasing temperature and vanishes above Tg. The poten- 
tial becomes a Debye screened Coulomb potential in the high temperature phase. 
While the leading high temperature behaviour as well as the structure of the heavy 
quark potential for temperatures well above Tg can be understood in terms of high 
temperature perturbation theory, it also is expected that non-perturbative effects 
like the generation of a magnetic mass term, cx g^(T)T (g(T) is a temperature 
dependent running coupling "constant"), in the gluon propagator may influence the 
spectrum in the deconfined phase. These non-perturbative effects in the magnetic 
sector will manifest themselves in correlation functions for the spatial components 
of gauge fields. 

(3-|-l)-dimensional renormalizable quantum field theories at high temperature, 
through dimensional reduction, can be reformulated as effective 3-dimensional theo- 
ries, with the scale of the dimensionful couplings given in terms of the temperature.^ 
In the case of SU(N) the effective theory includes interactions with adjoint matter 
(Higgs) fields, emerging from the temporal components of the original gauge fields. 
In particular, for SU(2) the effective action reads'*'® 
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being the three dimensional Wilson action. 
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Di{U)A4in) = Ui{n)A4n + ei)U+{n)-A4n) , 

describing the coupling to the Higgs field, and Sy being local Higgs self-interactions 
which are quadratic and quartic to their lowest orders, and carry additional powers 
of (jf^(T). For the Higgs field representing the temporal gauge field components, we 
have chosen the convention A4 = i^j-i A\<Tj where <Tj denote the Pauli-matrices. 

Basic properties of the gauge invariant correlation functions for spatial compo- 
nents of the gauge fields — e.g. the spatial Wilson loops — can be understood in 
terms of this effective theory. For instance, as this effective theory is confining, it is 
natural to expect that spatial Wilson loops obey an area law behaviour in the high 
temperature phase 

W{R,S) = {V[e^^^s^''^^^]) ^ e-"'^'^ , (9) 

where the (ordered) path is over a rectangle of area R x S, and ag is called the 
spatial string tension, although one should stress that it is not related to properties 
of a physical potential in the (3-|-l)-dimensional theory. In the case of QCD the 
effective theory itself is quite complicated even at high temperatures, as the non- 
static modes do not decouple from the static sector.^ An analysis of the temperature 
dependence of the spatial string tension, thus, yields information on the importance 
of the non-static sector for long-distance properties of high temperature QCD. 

The existence of a non- vanishing spatial string tension, ag, in the high tem- 
perature phase of (3-|-l)-dimensional SU(N) lattice gauge theory can be proven 
rigorously at finite lattice spacing.® However, despite its basic relevance for a better 
understanding of the non-perturbative structure of non-abelian gauge theories at 
high temperature, little effort has been undertaken to arrive at a quantitative de- 
scription of the properties of the spatial string tension. The first numerical study'* 
suggested that ag stays non-zero, but temperature independent, in the high tem- 
perature phase of QCD. Some indications for an increase of ag with temperature 
have been found recently.^ 

4. Measurements and Analysis 

We present here the results of a detailed, high statistics analysis of the spatial 
string tension. Some of our results have previously been published in Ref. 10. The 
finite (and zero) temperature SU (2) gauge theory has been simulated on lattices of 
size 32^ X Nr, with Nr ranging from 2 to 32 (Tab. 1). The simulations have been 
performed at two values of the gauge coupling, /3 = 2.5115 and /3 = 2.74, which 
correspond to the critical couplings for the deconfinement transition on lattices 
with temporal extent Nr^c = 8 and Nr^c = 16, respectively.® The lattice spacing 
thus changes by a factor 2.00 ± 0.04, where the error is caused by the uncertainty 
in both of the critical couplings. 

In each case we have performed between 2000 and 5000 thermalization sweeps. 
Measurements were taken every 100 sweeps. The integrated autocorrelation times 
have been calculated on some representative operators. In all cases they turned 
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Table 1: The simulated lattice volumes 32^ x Nr- The fourth column gives the 
number of gauge field configurations used in the analysis. Smtn labels the local 
potential Vt,s^j„ , used as an approximant for the asymptotic potential Vr (Eq. 10). 



out to be small compared to the above measurement rate. Besides the original 
data sets, we produced other sets by binning the data into blocks of two, and five 
measurements. The whole analysis was performed on the original data set as well 
as on the binned data. No systematic increase of errors with the block length 
was observed. We conclude that successive measurements can be treated as being 
effectively independent. 

Smeared on- and off-axis Wilson loops with the quark sources separated by 
multiples of the vectors (1,0,0), (1,1,0), and (2, 1,0) have been measured for the 
determination of the spatial string tension (correlator along z-axis). In addition to 
these separations we have choosen multiples of (1, 1, 1), (2, 1, 1), and (2, 2, 1) for the 
zero temperature cases. In each direction the quarks have been separated by up to 
half of the lattice extent which gives distances up to -j^N„ i^N„) in the finite 
(zero) temperature case. We apply the smearing technique presented in Ref. 11 
to the links orthogonal to the direction in which the correlator is taken (z for the 
spatial loops, and t for the "temporal" ones). 

After 40 such smearing iterations (temperature dependent) pseudo-potentials 
are determined from the Wilson loops VF(R, S) of shape R x 5, with R pointing 
into one of the above directions: 

y.(R) = £m^y.,.(R) , y-M^) = '- J^^s'}l) • (1°) 

Though, due to our smearing process, also temporal links are used for construction 
of the smeared Wilson loops, the (spatial) pseudo-potentials are not affected. This 
is due to the fact that the "3D-smeared" gg-creation operator as well as the corre- 
sponding "2D-smeared" operator share a nonvanishing overlap with all eigenstates 
of the Aig representation of the lattice symmetry group D^^h- Numerically, this has 
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Figure 1: The quality Q of various fits with different fit ranges [Rmm, Rmax\ against 
the fitted string tension k with error for T = Q, [3 = 2.74. In (a) all fits with 
confidence level a > 0.1 are plotted. In (b) fits with quality Q < Qmaxf^ are 
omitted. Note the different ordinate-scale. 

been checked for Nr = 2 at /3 = 2.5115. Because of faster convergence to maximal 
ground state overlap we have applied the 3D-smearing algorithm throughout our 
investigation. 

We use VT,s^jyj(R-) with finite Smtn as an approximant for the asymptotic 
value Vt(R-)^- Our requirement is that the systematic error caused by this cutoff is 
small compared to the statistical error on Vr • To guarantee this we choose Smtn in 
the way that for all S > Smtn the individual potential values as well as the fitted 
value for the string tension remain constant within errors. The latter criterion is 
important because a small but systematic shift of potential values may have an 
enhanced infiuence on certain fit parameters. 

The actual size of Smtn depends on statistics and ground state overlaps, as 
well as on the gap between ground state potential Vr and the first excitation 
in lattice units. Thus, for /3 = 2.5115 the asymptotics is reached earlier (in lattice 
units) than for /3 = 2.74. Also, for high temperatures not only the potential values 
are larger than in the zero temperature case but also the V^(R) — Vt(R-) gaps seem 
to increase. So, Smtn decreases with increasing temperature, too, as can be seen 
from Tab. 1. 

All our errors have been obtained from the scatter between (two times the num- 
ber of measurements) bootstrap samples. The whole analysis has been performed 

^ We prefer using "local masses" instead of exponential fits to Wilson loops with S > S^^j^ because 
of simplicity in data analysis. Due to the exponential increase of relative errors on pure gauge 
quantities with distance S, a fit would be governed by the values corresponding to S^^^ and 
Srrnn + 1> anyway. 
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Figure 2: Comparison between deviations of the potential at /3 = 2.74, T = 8Tc from 
the continuum symmetry parametrization V(R) = Vo + kR — c/R, and — Gl(R) 
(Eq. 12), which has been used to parametrize lattice artefacts. 



within each bootstrap sample, separately. For the zero temperature cases on each 
bootstrap a different covariance matrix between potential values Vt(R-) has been 
calculated (with a "subbootstrap"). A comparison with the result obtained by us- 
ing one and the same covariance matrix (obtained on the original data set) on all 
bootstrap channels did not show any significant change on errors. So, for the other 
fits the latter method was used for simplicity. All fits are correlated fits using this 
covariance matrix. We also cross-checked the results by comparing them with un- 
correlated fits since correlated fits may lead to totally wrong results if the statistics 
is too small to fix the covariance matrix accurately, or the assumed parametrization 
does not completely match the data.^^ 

The whole analysis has not only been performed for different Smtn on our 
original data sample and the two binned samples, but also with different fit ranges 
[Rrmn, Rmax]- For all these fits (about 1000) we have calculated the confidence 
levels a = 1 - (2'^°f^^T(dof /2))-^ dx x'^^^^^-^e-''^^ and defined qualities Q = 
a , — where k are the fitted values for the string tension. The fit with highest 
quality gives the "best" fit range. In Fig. 1(a) a histogram of the quality Q against 
the corresponding string tension value (with statistical error width) is plotted for 
all fits with a > 0.1 for the zero temperature, /3 = 2.74 case. In Fig. 1(b) the same 
is plotted, but only for qualities Q > Qmaxf^- We include the systematic error 
from the width of all the latter fits into all errors we state. As a check we also have 
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Figure 3: The pseudo-potentials Vt(R) minus the (constant) self-energy contribu- 
tions Vo (Eq. 11) on lattices of size 32^ x Nj- for /3 = 2.74 as a function of the spatial 
separation R, measured in lattice units. 
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Table 2: Summary of results from fits to the effective potentials using Eq. 11 on 
lattices of size 32^ x Nr ■ The values of Nr which correspond to the temperatures 
T/Tc are given in Tab. 1. In the case of the string tension the first error is purely 
statistical. All other errors include systematic effects due to the choice of the fit 
range in the way described in the text. 



defined a quality which uses the relative error of the Coulomb coefficient e instead 
of the string tension k. The resulting optimal fit ranges, as well as the systematic 
errors, have been almost the same in all cases. 

5. Results and Discussion 

5.1. Evaluation of the spatial string tension 

In the way described above, several fits to various parametrizations have been per- 
formed. For a three parameter fit with self energy Vq, Coulomb coefficient e and 
string tension k the values were found to be acceptable for Rmm > V5. Fits 
turned out to be more stable (and Rmm > V^) by use of the parametrization^^' ^'^ 



being Air times the lattice one-gluon-exchange for the infinite volume case. This 
last term takes account of the lattice artefacts present at small distances. The 
quality of this parametrization can be seen from Fig. 2 where a comparison (at 
T = 8Tc) between l/R — Gl(R), calculated from the potential data by use of the 
fitted parameter values and the theoretical l/R — Gl(R) is displayed. 

We also have experimented with lattice propagators on smaller and anisotropic 
volumes, replacing the integral by a discrete sum over the allowed lattice momenta 
but, in general, the quality of fits decreased with smaller volume. We conclude that 
the potential values seem to be insensitive to the low momentum cutoff. Reasons 
for this might be that the Coulomb interaction from the mirror sources is screened 
by the large (non-perturbative) linear part of the potential and that Wilson loops 
have a well defined inside and outside. 

Various other fit functions have been tried, including fits where the Coulomb 
part has been replaced by a logarithmic term, which would be expected in the three 
dimensional effective theory. We found that the fit parameter k = ago"^ only weakly 
depends on the actual parametrization of the short distance part of Vt(R)- Our 
results on the four fit parameters (Eq. 11) are summarized in Tab. 2. In Fig. 3 we 
display the potentials Vt(R), measured at different temperatures, for the /3 = 2.74 
case, together with the corresponding fit curves. We determine the spatial string 
tension in units of the critical temperature, ^/o^/Tc = \/kNj c, where Nj^^ = 8 (16) 
for [3 = 2.5115 (2.74). These numbers are given in the last column of Tab. 2. 

5.2. Scaling behaviour 

A calculation of the "zero" temperature string tension in the confining phase on a 
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Figure 4: The physical zero temperature potentials minus the self energy parts Vq 
(Eq. 11) at /3 = 2.74, and /3 = 2.5115 in units of the /3 = 2.74 lattice spacing. The 
l3 = 2.5115 values have been scaled by a factor 0^=2.5115/0/3=2.74 = 2.016. 



lattice of size 16 x 32^ (324) at /3 = 2.5115 (2.74) yields 

r- _ r 0.1836± 0.0013 ,/3 = 2.5115 . . 

^"^^ ~ \ 0.0911 ± 0.0008 ,/3 = 2.74 ' ^ ' 

which corresponds to a change in lattice spacing 0^=2.5115/0/3=2.74 = 2.016 ± 0.023, 
and is consistent with the factor 2.00 ± 0.04, obtained from the calculation of the 
critical couplings for the deconfinement transition.® This is visualized in Fig. 4 where 
the zero temperature potentials, measured at the two different /3 values, are scaled 
by this factor 2.016. 

In Fig. 5 we compare the spatial string tension calculated at /3 = 2.5115 and 
2.74 at different temperatures. We find that — as in the zero temperature case — 
our data sets are consistent with each other down to a temporal extent as small 
as Nr = 2. Thus, scaling violations in the ratio ^/o^/Tc are negligible. This 
demonstrates that the spatial string tension, indeed, is relevant to high temperature 
QCD as it persists in the continuum limit. Moreover, ag coincides with the physical, 
zero temperature string tension for T < Tg. 

5.3. Comparison with the effective theory 

The coupling of the Yang-Mills part of the action of the effective 3-dimensional 
theory (Eq. 6), §3, derived from a (3+l)-dimensional SU(N) gauge theory at high 
temperature, is given in terms of the temperature and the four-dimensional coupling 
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Figure 5: Square root of the spatial string tension in units of the critical temper- 
ature versus temperature, calculated at two different values of the gauge coupling. 
The broken line gives the result for the ratio of the physical string tension to the 
deconfinement temperature averaged over several values of the critical coupling.® 
The horizontal error bars indicate the uncertainty in the temperature scale Tg. 

g{T) as (jf| = g^(T)T. Although the temperature will also set the scale for the 
Higgs self-couplings in the 3-dimensional theory, these couplings will in general 
have a different dependence on the four-dimensional gauge coupling g(T)}^ The 
functional dependence of <Js{T) on g^(T) and T, thus, is not apparent from the 
general structure of the effective action. 

However, in a pure three-dimensional SU(N) gauge theory dimensionful quan- 
tities are proportional to an appropriate power of the three-dimensional coupling 
gs. If the temperature dependence of the pure gauge part of the effective action 
dominates the spatial string tension we would expect to find 

^a,(T) = cgl = cg^(T)T , (14) 

where the temperature dependent running coupling constant g^(T) should, at high 
temperatures, be determined by the perturbative /3-function of SU(N) in four di- 
mensions. 

In Fig. 6 we have plotted T/ ^(7s{T) against T. From Eq. 14 this ratio is 
expected to be proportional to g~^(T). We have fitted these data to the two-loop 
formula for the coupling in SU(2) gauge theory with the scale parameter At, 

9-'(T) = T^lnT/A^ + -^ln(21nT/AT) . (15) 
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Figure 6: The ratio of the critical temperature to the square root of the spatial 
string tension versus temperature for /3 = 2.74. The line shows a fit to the data in 
the region 2 < T/Tc < 8 using the 2-loop relation for g(T) given in Eq. 15. 



We find that the temperature dependence of the spatial string tension is well de- 
scribed by Eq. 14 and Eq. 15 for temperatures above 2Tc. From the two parameter 
fit to the data shown in Fig. 6, we obtain 

V<T.(T) = (0.369 ± O.OU)g''iT)T , (16) 

with At = 0.076(13)Tc. We note that the second term in Eq. 15 varies only little 
with temperature. A fit with the one-loop formula thus works almost equally well; 
it yields a!^^^ = 0.050(10)Tc and c = 0.334(14) for the coefficient in Eq. 16. 

It is rather remarkable that the spatial string tension depends in this simple 
form on the perturbative SU(2) /3-function already for T > 2Tc. Moreover, we find 
that even quantitatively the spatial string tension agrees well with the string tension 
of the three-dimensional pure SU{2) gauge theory, -'^'^ = (0.3340 ± 0.0025)(jf|.+ 
We take this as an indication that, indeed, the spatial string tension is dominated 
by the pure gauge part of the effective three-dimensional theory. There are also 
numerical indications^^ that a Higgs coupling of the strength suggested by per- 
turbative arguments® has no significant infiuence on the potential values of three 
dimensional Yang-Mills theory. We note that the value for g^(T), determined here 
from long distance properties of the (3-|-l)-dimensional theory, is about a factor two 

+ We have independently confirmed tliis value by use of smeared Wilson loops instead of Polyakov 
line correlators.-^^ 



larger than what has been obtained by comparing the short distance part of the 
(3+l)-dimensional heavy quark potential with perturbation theory.® 

Another quantity characterizing the high temperature behaviour is the electric 
Debye screening mass, rUei. It has been computed to one loop,"'^® rUei = \/2/'ig{T)T , 
with g{T) being a running coupling constant with scale parameter Kej. In one ap- 
proximation this scale parameter has been computed to be"'^^ Kgi = Q.Qb^Kjj^, which 
can be compared to our result from the spatial string tension, -* = {].{](S2{lA)Kjj^ 
(here we have used® Tg/Aj^ = 1.23 ± 0.11). It should be emphasized, though, 
that there is no unique definition of a temperature dependent coupling, g(T), but 
it is tantalizing that g(T) given by Eq. 16 might be quite close to g(T) determined 
from the Debye screening mass. We intend to investigate this possibility in the near 
future. 
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